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We introduce the simplest model to describe parametric interactions in a quadratically nonlinear 
optical medium with the fundamental harmonic containing two components with (slightly) different 
carrier frequencies [which is a direct analog of wavelength-division multiplexed (WDM) models, 
well known in media with cubic nonlinearity] . The model takes a closed form with three different 
second-harmonic components, and it is formulated in the spatial domain. We demonstrate that the 
model supports both polychromatic solitons (PCSs), with all the components present in them, and 
two types of mutually orthogonal simple solitons, both types being stable in a broad parametric 
£\j | region. An essential peculiarity of PCS is that its power is much smaller than that of a simple 

■ (usual) soliton (taken at the same values of control parameters), which may be an advantage for 

experimental generation of PCSs. Collisions between the orthogonal simple solitons are simulated 
£Sj ' in detail, leading to the conclusion that the collisions are strongly inelastic, converting the simple 

solitons into polychromatic ones, and generating one or two additional PCSs. A collision velocity at 
which the inelastic effects are strongest is identified, and it is demonstrated that the collision may 
be used as a basis to design a simple all-optical XOR logic gate. 
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Wave mixing at different carrier frequencies, of which generation of higher harmonics is a well-known example, has 
had a long history of investigation g|. In an optical medium whose symmetry group lacks a center of inversion, the 
r— h ■ lowest-order nonlinear response is quadratic, which gives rise to three-photon interactions. The resultant three-wave 
i " mixing is resonant when the condition UI3 = oji + 0J2, imposed on frequencies of the three waves, is met. In the 
particular case when u>i = LJ2, the process reduces to the second harmonic generation (SHG). 
' Solitons in quadratic [x^] media has been a major area of research recently (see Ref. B for a review). The first 
experimental observation of x^ solitons was reported for type-I SHG, which involves exactly one component of each 
harmonic, in the (2 + l)-dimensional (bulk) geometry Q. The observation of solitons in the (1 + l)-dimensional 
geometry (planar nonlinear waveguide) followed shortly afterwards Q . Much theoretical work has been performed for 
the solitons in both type-I and type-II SHG, the latter case involving two components of the fundamental harmonic, 
corresponding to different polarizations, and a sing le component of the second harmonic §, @, g, g The SHG 

£Nl ■ process in an isotropic medium, with two polarizations at both harmonics, has been considered too [|12[ 

More complex cases of multi-resonance wave mixing in quadratically nonlinear media have not received much 
attention because of serious difficulties with their experimental realization using the birefringence-based wave-vector 
matching schemes. However, the recent rapid development of the quasi-phase-matching (QPM) technique has changed 
the situation. The technique, originally proposed long ago ||, relied upon periodic structures (usually, periodically 
poled ones) with an alternating sign of the quadratic nonlinearity. Recently, the QPM technique has been extended 
. ! from periodic to Fibonacci-series-based structures jl3| , and further to fully quasiperiodic ones jL4| (see also Ref. jl5| ) . 
This makes it possible to essentially relax conditions on the wave-vector difference for the matching to take place. 
Results presented in Ref. |0 show, both theoretically and experimentally, not only high values of effective x^ 
coefficients, but also that one can attain the wave-vector matching simultaneously for several sets of waves involved 
in the nonlinear interactions. In particular, the development of the QPM technique has been an incentive to study 
multi-resonance systems jl^, [H], |l8| |l9|], a possible application of which may be design of soliton-based logic gates. 
Indeed, a soliton may naturally be used as a bit of information, and the interactions of solitons can potentially support 
logic operations. 

In this work, we develop a model of nonlinear mixing between two fundamental-harmonic waves with different 
frequencies in a quadratic medium. Via the x nonlinearity, they generate three different wave components of the 
second harmonic. Note that interactions between waves with different frequencies in optical media with cubic (x ) 
nonlinearity is a well-known topic, which has extremely important applications to the WDM multichannel format 
of data transmission in fiber-optic telecommunication links p0[ |; however, it appears that the issue has not yet been 
studied for x^ media, in which the mixing may be realized in both spatial and temporal domains. 



We formulate the model in section II, and produce its stationary soliton solutions in section III. These may be both 
fully polychromatic solitons (PCSs), including all the five field components, and special (simple) solutions of "A" and 
"B" types, which amount to ordinary two-component SHG solitons in mutually orthogonal (non-intersecting) subsets 



of the five- wave system. In section IV, wc demonstrate, by means of direct simulations, that the solitons of all these 
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types are stable in a broad region in the system's parameter space. In the same section, we simulate collisions between 
the simple solitons. When they overlap, the nonlinear interaction generates a component which was absent in both 
of them before the collision, which makes the final result of the collision strongly inelastic: the former simple solitons 
develop extra components and become polychromatic ones. Additionally, one or two extra solitons are generated by 
the collision. We demonstrate that the collision between the initial simple solitons may be a basis for an all-optical 
XOR logic gate. The paper is concluded by section fv|. 



II. THE MODEL 



We consider the interaction of five waves in a diffractive dielectric medium with the quadratic nonlinear suscep- 
tibility. The carrier frequencies of the waves satisfy resonant conditions, uj 3 = 2u>x, lo 4 = lui + uj 2 , = 2^2, so 
that uj\ and uj 2 may be classified as carriers of the fundamental-harmonic components, while u> 3 , L04 and W5 represent 
three components of the second-harmonic wave group, generated by the two fundamental-harmonic components via 
the nonlinearity. 

Assuming, as usual, that the wave envelopes Ex,E% and E 3 ,E 4 , E§ of these components are slowly varying ones, a 
system of five equations, coupled parametrically through the components of the x^ nonlinear susceptibility tensor, 
can be derived from the Maxwell's equations to govern the evolution of the waves in the spatial domain (the derivation 
follows the well-known procedure worked out for the usual type-I and type-II x^ interactions, see a detailed account 
given in the review 

pij? f) 2 K 

2ik i-iT- + "inr + 2X1^* i^e^ 1 ' 1 + V2x 2 E* 2 E 4 e* Ak ^ = 0, 
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where k\ through k§ are the corresponding carrier wave numbers, z and x are the propagation and transverse coor- 
dinates, Akimn = fa — k m — k n is the wave- vector mismatch, and Xn = 87rw^/c 2 x' 2 \ x^ being an element of the 
quadratic susceptibility tensor. The extra factors of 2 in front of the coefficients xi an d X.3 m the first two equations 
reflects the fact the equations may be derived from a Lagrangian, and the factor %/2 in front of all the terms containing 
Xi is added by definition, to simplify subsequent rescalings. 

Equations (fy) are normalized by measuring x and z, respectively, in units of the input beam size tq and diffraction 
length Zd — r^k^ at the frequency lo\ + u>2- Introducing dimcnsionless fields 

Ex = (U1/V27) exp (ifcz) , E 2 = (u 2 /V2 7 ) exp [ifoz) , E 3 = (1*3/27) exp (2ifaz - iAk 311 z) , (2) 
E 4 = (U4/V27) exp (*(/3i + fo)z - iAk 412 z) , E 5 = (u 5 /2 7 ) exp (2i(3 2 z - iAk 522 ) , (3) 

and 7 = XiTq, a normalized system of equations is obtained: 

h ( ,dui \ d 2 Ul 

kl \ ~b~z ~ 1 J ~dx 2 XlUlU3 + ^ 2U 2«4 = 0, 

k 2 ( .du 2 \ d 2 u 2 

k A \ ~dz~ ~ J ~dx T X2UlU4 + ^ 3M 2«5 = 0, 

ki ( .8114 . . \ d 2 u A 

2-j^ y-fcT - \"l +P2 - Afc4i 2 )U4 I + -^j" +X2U1U2 = 0, 
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where Xn = Xn/xi ( n — 1) 2, 3, i.e., X i = 1)) fti,2 being two phase-velocity shifts. To reduce the number of parameters 
in Eqs. (|J), the fields and coordinates can be rescaled further. Defining u n — PiU n , z — Z/ftx , and x = X/ \/|/3i |, 
we obtain 



arr FPU 

i-q7T + -q^-u 1 U 1 +xiUZU 3 + X2U;U 4 = 0, 



dU d 2 U 

i-Q^ + —^-a 2 U 2 +X2UZU 4 + x 3 U;U 5 = 0, 



qjj d 2 U 
2i ~di + ~8X^ ~a 4 U 4 + X 2U 1 U 2 = 0, 
.dU, d 2 U, „ 1 



2i--g + ^-a 5 U 5 + - Xs Ul = 0. 



where a x = 1, a 2 = fo/fa, «3 = (4/3i - 2Afc 3 n)/^i, «4 = (2/?i + 2/3 2 - 2Afc 4 i 2 )//3i, and a 5 = (4/3 2 - 2Afc 522 )/A- 
We assume that |/c4 — fc 2 — fci | <C k n , and everywhere, except for the phase- mismatch parameters introduced above, 
one may set fc 3 = k 4 = k 5 , i.e., the three components of the second- harmonic wave have similar wave numbers. 

Equations (|5|) assume that the three parametric interactions ("vertices"), which couple, respectively, the wave 
triplets (1,1,3), (2,2,5), and (1,2,4), are nearly phase matched (as it was mentioned above, a possible way of achieving 
this may be provided by the QPM technique based on quasiperiodic structures [Q). It is straightforward to see that, 
like the model describing the type-II SHG (|, ||, |l0|], Eqs. (||) have two Manley-Rowe invariants, namely, the total 
power, 



/oo 
|C/i| 2 + \U 2 \ 2 + 4 (\U 3 \ 2 + \U 4 \ 2 + \U 5 \ 2 ) dx, (6) 
-oo 

and the power imbalance, 

/oo 
|E/i| 2 - |[/ 2 | 2 +4(|[/ 3 | 2 - \U 5 \ 2 )dx. (7) 
-oo 

The present model does not include walkoff terms (group- velocity mismatch). While a detailed discussion of the 
walkoff is beyond the scope of this work, it is relevant to mention that QPM and similar techniques, such as tandem 
structures pi] , make it possible to suppress the walkoff ^2j. We also note that QPM can give rise to an effective 
cubic nonlincarity psjj , which may compete with the underlying quadratic interactions []l5| , |2^| . For this reason, cubic 
terms should sometimes be added to a dynamical model, but this is not an issue for immediate consideration in the 
present context. As concerns the physical realization of the system, estimates using typical values of the relevant 
physical parameters in such quadratically nonlinear materials as lithium niobate and KTP |l4j , p"5| ] suggest that a 
necessary (quasi) period of the QPM structure is ~ 10 /im, and the power and transverse size of the soliton beam 
are expected to be ~ 20 /jm and 100 mW, respectively. 

If still more resonances are allowed, other essential wave components may appear, for instance, those corresponding 
to the combinational frequencies lui 2 = 2lo\ — lu 2 and lo 2 i = 2lo 2 — io\ (obviously, they belong to the fundamental- 
harmonic wave set). Denoting the corresponding wave numbers as k\ 2 and k 2 \, we see that these new components will 
indeed be essential if the wave-number triplets (12,2,3) and/or (21,1,5) are nearly phase-matched. The accordingly 
modified system will include seven components (four pertaining to the fundamental harmonic, and three to the second 
harmonic) and five vertices. However, such a situation seems much more exotic (five simultaneous resonances) than 
the more generic possibility of three simultaneous resonances underlying the model considered in the present work. 



III. STATIONARY SOLITON SOLUTIONS 



Particular exact solutions of Eqs. (|5|) for stationary solitons can be sought for in the form 

U n = ^„sech 2 (XX) ,n=l, 5, (8) 

where A n are amplitudes, and A is the inverse width of the soliton. Inserting this (|§|) into Eqs. (|^), a solution can be 
obtained if a± = a 2 — a 3 — a 4 = a$ = a and A = i/a/2. The amplitudes are found to be 

fiA 2 v A 2 
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A 4 = M = ^, (9) 

where ^(xi/2-x|)/(x!/2-x|). 

In addition to the exact solution based on Eqs. (g) and (Q), other particular solutions can be found setting t/2,4.5 = 
or, alternatively, C/1,3,4 = 0. Equations (||) then reduce to the well-known type-I SHG model ||, [| fr), 

i_ + ^-V + V W = 0, (10) 
dW d 2 W 1 n 

where V and are the fields at the fundamental and second harmonics. In this case, there is a well-known special 
exact solution corresponding to p = 1 in Eq. (O) [E5j , 



V = (3/V2) scch 2 (X/2) , VF = (3/2) sech 2 (X/2) . (12) 

When p / 1, a family of stationary soliton solutions to Eqs. jl0| ) and ( pi]) can be found numerically [Q (or approxi- 
mated analytically by means of the variational method |2q| ). We will refer to the general solution for the case when 
Ui = V, U3 = W, and E/2,4,5 = as an "A" type soliton, while the opposite "B" type is defined as the one with 
U2 = V, U5 = W, and J/1,3,4 = 0. Both these types will also be called "simple" solitons. 

General solutions for the polychromatic (five-wave) soliton can be found from the stationary version (the one with 
dU n /dZ = 0) of Eqs. (^|) by means of the standard numerical methods for two-point boundary- value problems. In 
Figs, [fl a nd 0, comparison is made between the general five-wave solitons and the particular solutions generated by 
Eqs. po| ) and ([ll]) at equal values of all the parameters. In Fig. [j], one can see that the A soliton is much larger in 
amplitude, while its width is not widely different from that of the polychromatic one. As a consequence of this, the 
power Q [see Eq. (^J)] of the A-type soliton shown in Fig. [l is Q = 36, while for PCS in the same figure, Q = 8 [the 
other invariant is Qi m b = in the case considered, see Eq. (7|)]. 

This drastic difference in the powers can be understood: in the case of the full PCS, one has two nonlinear terms in 
the first two equations (^), rather than one term in the case of the simple solitons; therefore, the amplitude necessary 
to compensate the spreading out of the beam due to the diffraction term is, roughly, twice as small in comparison 
with the simple solitons, or, eventually, the power is ~ 4 times as small. To check whether the power of the PCS is 
indeed essentially lower than that of simple solitons in the general case, in Fig. |^ we show a typical example of the 
change of the two powers with the variation of a.2 or, equivalently, p in Eq. (|l l|) . As is seen, PCS indeed persistently 
maintains a lower power Q than its A-type counterpart. On varying the different parameters a n ,Xn the value of 
Q of the polychromatic soliton may be increased (or decreased) in value from that portrayed in Fig. || but it only 
exceeds that of the simple solitons when \2 < 0.7. In the limit of \2 — ► Eqs. (||) decouple and the polychromatic 
soliton tends to a "double" simple soliton i.e. both "A" and "B" type solitons existing in one envelope. The double 
simple soliton by definition has twice the power of single simple soliton. Provided the nonlinearity coefficient \2 is 
large enough (x2 > 0.7) it is possible to conclude that PCSs in a quadratically nonlinear medium may be produced 
from waves with different frequencies at a much lower net input power than the ordinary SHG solitons, i.e., it may be 
essentially easier to generate PCSs in the experiment and use them in potential applications. Of course, these results 
are meaningful provided that these solitons are stable. 



IV. STABILITY AND INTERACTIONS OF THE SOLITONS 



A. The stability of the polychromatic solitons 

To test the stability of PCSs, we solved the full system of Eqs. (||), employing the numerical split-step fast-Fourier- 
transform method. The simulations were performed with a computational grid of 2048 points, the transverse and 
propagation step-sizes being, respectively, h x 0.02 and h z = 0.01. The integration domain had the transverse size 
100 (—50 < x < +50), which is by far larger than any x size relevant to the solitons, see Figs. 4, 5, and 9 below. 
Absorbing layers were placed at the edges of the computational domain to prevent reentering of radiation. Specially 
monitoring interaction of the radiation waves with the absorbers, we have verified that, in all the cases considered, 
no reflection took place indeed. 

Evolution of initial configurations close to the stationary solutions was simulated for a variety of parameters. To 
impose initial perturbations, values of the initial amplitude and width of the pulses were altered against the exact 
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stationary solutions. From the results of the simulations, we have concluded that PCSs survive, clearly remaining 
stable, as long as the simulations were run, the maximum simulation length being 300 diffraction lengths of the soliton. 

To further test robustness of PCSs, we ran numerical experiments in which the solitons were successfully generated 
from initial Gaussian pulses, launched in the pump fields U± and U2 (with no initial second-harmonic components, 
which corresponds to the generation of SHG solitons in real experiments ^7j). The results demonstrate that PCSs 
not only are stable against small perturbations, but also play the role of strong attractors in the system. An example 
of PCS generation from the Gaussian beams is illustrated by Fig. |3|. Typically, there is a period of strong relaxation 
of the beams, where the amplitudes fluctuate and excess energy is radiated away. The second-harmonic components 
^3,4,5 are generated, and the soliton arranges itself to a quasi-stationary form, which then propagates, in a stable 
fashion, over a distance in excess of 100 diffraction lengths (as long as the simulations were run). As a result of many 
runs of systematic simulations, PCSs have been found to be attractors in a broad range in the system's parameter 
space. We concentrated mainly on the parameter space a n = 1..5 and \n = 1-2 where the PCS definitely had a lower 
value Q than the simple solitons. In connection to this, it is relevant to note that the above-mentioned families of 
the simple A and B solitons are also stable in the ordinary SHG model, except for a small region of their existence 
domain J|, 0] . 

B. Collisions between orthogonal simple solitons 

An interesting possibility is to consider collisions of the mutually orthogonal simple solitons of the above-mentioned 
A and B types. The overlapping between the colliding solitons will give rise to the generation of the field E/4, which 
is absent in both A and B solitons in their pure form, and the issue is how the generation of this field will affect the 
interaction between the solitons. Equations ( |TT| ) have the property of the Galilean invariance, so "moving" solitons (in 
fact, the solitary beams propagating at an angle in the planar waveguide) can be constructed by the transformation 

V(X,Z) = V Q (X -CZ)exp(iCX/2~iC 2 Z/4) , (13) 
W{X,Z) = W Q (X -CZ)cxp(iCX -iC 2 Z/2) , 

where C is the "velocity" (actually, a slope) of the "moving" soliton. 

We collided the A and B solitons with an initial separation between their centers Xq = 20, varying their velocities 
±C. A representative set of values of the parameters for which the results are displayed below are ai .,5 = 1, 
Xi = X3 — 1; an d X2 = 2; many simulations run for other values gave quite similar results. 

In the case of moderate collision velocities (and zero phase difference between the solitons) , the generation of the 
field Ui in the course of the collision gives rise to a third polychromatic soliton with the zero velocity. Trajectories 
of the initial A and B solitons alter, and they appear from the collisions as PCSs too, i.e., the collision adds to them 
components which were initially absent. The three post-collision solitons have roughly equal powers, with a mirror 
symmetry in the power distribution amongst individual components, see Fig. ||. At higher collision velocities (again, 
for the zero phase difference), four PCSs are generated by the collision, all having a nonzero velocity (see Fig. |[). As 
the collision velocity increases, less and less interaction takes place, until the initial solitons pass through one another 
unchanged (elastically) . In all these cases, 90% of the initial power is typically converted into the resulting PCSs, i.e., 
radiation loss is not a dominant actor in the collision dynamics. 

Results produced by many runs of the simulations for the collision of the in-phase (zero-phase-difference) solitons 
are summarized in Figs. ^ and 0. The constant initial separation xq between the solitons in the simulations means 
that the increasing velocity corresponds to an increasing incidence angle 0\. As may be naturally expected, the 
collision alters the trajectories of the solitons most in the case of low velocities. At higher velocities, the outer solitons 
keep essentially the same velocity after the collision as they had before it, while the two additionally generated inner 
solitons are significantly slower. 

The most significant inelastic effects occur when the solitons collide with the velocities ±0.4. In this case, the exit 
trajectory is altered the most (see Fig. |6|), and, as per Fig. fj], the largest part of the net power is transferred into the 
newly generated harmonic components of the outer solitons. 

The interaction between the A and B solitons produces nontrivial results also in the case when they have zero initial 
velocities, but their tails overlap at the initial position. For the initial separation xq — 5, the result of the interaction 
with C = is displayed in Fig. ||, in the form of the distribution of the fields U\ and U2 produced by the interaction. 
Again, three PCSs appear, the central one with the zero velocity, and two outer solitons with nonzero velocities. In 
the course of the evolution, a phase difference develops across the component fields. In fact, it is this phase difference 
which repulses the outer solitons, lending them a nonzero velocity. 

It is well known from theoretical |^8| and experimental [^7| studies of collisions between SHG solitons of the usual 
type that the result strongly depends on their relative phase at the collision point (provided that the colliding solitons 
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have nearly identical amplitudes): they attract each other and therefore interact strongly in the case of the zero 
phase difference, and repel each other if the phase difference between the fundamental-harmonic components is n. 
However, in the present system the effect of the phase difference on collisions between the simple solitons of the 
orthogonal types is much weaker. Simulations performed with various phase differences (including it/2) between 
the two orthogonal simple solitons does not display any change in the post-collision dynamics of the beams when 
compared to the zero phase difference case. Also it may be seen from the structure of the underlying system (|^) 
that phase differences between the "A" and "B" solitons will have little effect. Indeed, these equations are exactly 
invariant against the transformation u\ — > iti exp (i<^i), it 3 — > U3exp(2i0i), ui — > U2 exp (1^2), u$ — > U5 exp {2i<f>2) 
and M4 — > 114 exp (itfii + ifo) with an arbitrary phase shift 4>i,2- 

The interactions between the simple solitons in this system may find application as a basis for an all-optical logic 
gate. Indeed, consider two solitons of the "A" and/or "B" types, propagating in such a way that they will collide. 
If the solitons belong to the same type (i.e., the configuration is AA or BB), then they attract or repel one another, 
depending on their relative phase p8jj , or, if the relative velocity is high enough, they simply pass through one 
another. But if they belong to the opposite types (an AB configuration), at least three PCSs are produced by the 
collision even at high velocities. This is a behavior which is expected from an exclusive OR gate, alias XOR. The 
actual outcome can be easily established by checking the C/4 content in the output (by means of an appropriate optical 
filter). An advantage of such a design of the XOR gate is that any output beam will be a stationary soliton (even in 
the case of a strongly inelastic collision), which makes it convenient for further manipulations (cascadability) . 



V. CONCLUSION 



We have introduced the simplest five-component model of polychromatic solitons (PCSs) in a quadratically nonlinear 
optical medium. The existence and stability of both polychromatic and two types of simple solitons have been 
demonstrated. An essential peculiarity of PCSs is that their power is much smaller, at the same values of the 
control parameters, than the power of the usual two-component (simple) solitons. We have also performed systematic 
simulations of collisions between mutually orthogonal simple solitons, concluding that the collisions are strongly 
inelastic (including the interaction between two solitons with the zero initial velocity), giving rise to transformation 
of the simple solitons into polychromatic ones, and generation of one or two additional PCSs. A value of the relative 
velocity at which the inelastic effects are strongest has been found. We have also shown that the collision may serve 
as a basis to design a simple all-optical logic gate of the XOR type. 
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FIG. 1: Profiles of a polychromatic soliton (solid lines) and a type "A" soliton (dotted). Common values of the parameters for 
these solutions are ai,...,5 = 1, Xi = lj X2 — 2, and X3 — 1- 
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FIG. 2: The power invariant Q for both the "A"-type (upper line, vs. p) and full polychromatic (lower line, vs. a.2) solitons. 
Values of the parameters are the same as in Fig. El. 
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FIG. 3: The generation of a polychromatic soliton from the Gaussian input U\ = 2e~ x 1/4 and U2 = 3e~ x After an initial 
period of strong relaxation, a weakly oscillating soliton is produced. The evolution of the amplitudes of the U n components vs. 
Z is plotted. The parameters used in this simulation were ai = Q2 = 1, 03 = 04 = 05= 4, xi = X'i — 1) an d X2 = 2. The 
power invariants are Q = 13\/27r and Qimb = — 5\/2tv. 




FIG. 4: The collision between simple solitons of the mutually orthogonal A and B types, with the initial velocities ±0.3 and 
zero phase difference. The parameters are a n = 1, xi = Xs = 1> X2 = 2, and C = 0.3. The field components are displayed 
in the panels: Ui (a), U2 (b), U3 (c), U4 (d), and C/5 (e). The panel (f) shows a combined contour plot of the Ui and U2 
amplitudes. Note the mirror symmetry of the profiles of the components U\ and Uz, which is a consequence of the fact that the 
imbalance invariant Qimb is zero in this case. It is obvious that the simple solitons become polychromatic after the collision, 
and the third polychromatic soliton with the zero velocity is generated. 




FIG. 5: The same as in Fig. g| except that the collision velocities are larger, ±1.0. A zero- velocity soliton is no longer generated 
after the collision. Instead, four solitons are observed in the post-collision state. In the panel (f) the outer beams terminate 
prematurely because they hit the absorbing sponges used in the numerical scheme. 
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FIG. 6: The exit angle (in degrees) 8 C for the outer solitons, as found from the simulations of the collision, vs. the incidence 
angle 9i. The angles are between the solitons' trajectories and the z axis in the (x, z) plane, so that an increase in the incidence 
angle corresponds to an increase in the collision velocity. Naturally, the collision becomes less inelastic as 0i increases, and 
eventually the exit angle becomes nearly equal to the incidence one. 
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FIG. 7: Distribution of the powers Q n = \U n \ 2 dx between the components of the former A or B soliton after the collision, 
as a function of the soliton velocity before the collision. At higher velocities, the collision gives rise to 4 (rather than 3) solitons, 
which take a part of the energy. 
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FIG. 8: The combined amplitude contour plot of the pump waves U\ and U2 in the case when the initial velocities of the simple 
solitons "A" and "B" are zero. 



